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1. Introduction 

A number of questions raised during this meeting can be clarified by the 
study of the emission of SNe la at late times. Do WDs explode below the 
Chandrasekhar mass? Or, on the contrary, are the differences that we see 
among Type la supernovae (SNe la) just variations within the paradigm of 
Chandrasekhar-mass exploding WDs? The evolution of the spectra along 
the exponential decline of the SN light curve, and the light curve itself at 
those phases, bring us information on the density profiles, composition gra- 
dients, and mixing within the supernova ejecta. Thus, important clues on 
nucleosynthesis, hydrodynamics, and chemical stratification can be inferred 
from the study of the late emission. Indeed, the initial steps on this path 
have already probed the underlying explosion model for thermonuclear su- 
pernovae well beyond a mere rough assessment: the studies of Meyerott 
(1978, 1980) and Axelrod (1980) provided a firm support to the exploding 
CO WD model for SNe la. Late-time studies of the other type of super- 
novae, the core-collapse ones, have also been crucial for their theoretical 
understanding. Fransson & Chevalier (1989) discussed the evidence of mix- 
ing in Type II supernovae such as SN 1987A from its effect in the late 
emission. Their analysis of the emission in the explosion model of CO WDs 
accreting from He stars revealed, as well, that it was inadequate to account 
for Type lb supernovae. The model based on massive progenitors was, on 
the contrary, favored by their calculations. An essential reference to these 
and other issues on late-time emission is Fransson (1994). 

In the present review, we will concentrate on the physical processes 
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which play a role in the formation of the spectrum and light curves of SNe 
la at late phases, and on the utility of those phases to clarify the way in 
which C + WD explode. The role of e"*", 7-rays, the paths of the high- 
energy radiation to thermalization through a number of complex cascades, 
the final optical display and the information contained in it, makes super- 
novae at these phases a physical laboratory dealing with a large variety of 
radiative processes, with energies ranging from MeV down to a few orders 
of magnitude below one eV. We will show how it is possible to come out 
with first-order answers, in the way radiation escapes from the labyrinth 
of plasma processes. 

2. Late— time approach: how late? 

The mere understanding of the spectral evolution of a supernova is a goal 
by itself. The temporal development of the spectral features and a proper 
identification of some of them are still a pending task. We can benefit, 
however, of some simplification if we decide to study a particular phase 
of the spectral evolution. Approaches to study the photospheric phase are 
discussed in this volume by various authors. Here, we concentrate on the 
times when the continuum emission has dropped substantially and nebular 
emission increases. We will model the SN at those phases; later on they 
will become the subject of study of our SN remnant colleagues, as the 
density becomes significantly lower and the radiative shocks originated by 
the interaction of the ejecta with the interstellar medium play an important 
role. The nebular phase starts in SNe la at about 80 days after explosion. 
At this age the ejecta become optically thin and a photosphere does not 
longer exist. This transition proceeds in a continous way and the continuum 
emission drops while the the intensity of nebular emission from forbidden 
transitions of Fe-peak elements grows larger. From the condition that the 
optical depth to continuum processes becomes negligible, Kcpr ^ 1, and 
taking into account that the supernova expands homologously with r = 
f X it can be seen that 



where Kc is the continuum opacity, M is the total mass, and vg is the 
velocity in units of 10^ cm s Thus, about 250-300 days after explosion 
the spectrum is really nebular. 

At this phase the steady state hypothesis for ionization balance and 
thermal balance hold (see §6). This allows some simplification of the com- 
plex physics of nonthermal processes in these non-LTE plasmas. The optical 
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depth for absorption in permitted transitions of the emitted photons is also 
significantly larger: 

Tbf ^Gbf nr (2) 

and it is of the order of 

nj ^ 6.9 X 10^ auv 7^ ^ 1 X 10^ (3) 

(vgty (^9)2 

Bound-free transitions will be followed by recombination cascades down 
to the ground state. Ionization resulting from this recombination plays a 
significant role. 

3. The energy source: radioactivity from ^^Co 

The radioactive decay 

56Co Fe 

provides the long-lasting luminosity source at this stage. Such a decay gives 
rise 81% of the time to a 7-ray photon and 19% to a e"*". 7-ray photons are 
emitted with a spectrum of energies reaching up to 1.4 MeV. They bear 
about 96.5 % of the energy of the ^^Co decay. The emitted positrons have an 
energy spectrum extending up to the endpoint kinetic energy Emaa;= 1-459 
MeV, and they account for 3.5% of the energy of the ^^Co decay. 

Compton scattering of the emitted 7-rays gives rise to a pool of non- 
thermal electrons called primaries. The optical depth to 7-rays is given by 
r-y ~ GK^UeR, where (Ta'jv is the Klein-Nishina cross section for Compton 
scattering. Once 7-ray photons suffer Compton scattering either they do 
not lose any significant amount of energy (forward scattering), or they lose 
significantly their energy, becoming readily absorbed (unable to produce 
further energetic primaries). This has suggested (Sutherland & Wheeler 
1984) the adequacy of treating the Compton scattering process as an ab- 
sorption process, for applications related to the energy deposition of 7-rays. 
Swartz, Sutherland, & Harkness (1995) review the result of comparing the 
purely absorptive approach and the Monte Carlo one. The detailed geome- 
try of an arbitrary 3D distribution of radioactive material can be accounted 
for within the first approach, saving CPU time as compared with 3D Monte 
Carlo calculations. However, if one needs to calculate the energy spectrum 
of the 7-rays and that of the primaries created by Compton scattering, 
Monte Carlo calculations are advisable. 

The outcome of the calculations is the fraction of radioactive energy 
deposited in the supernova ejecta: 
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Figure 1. Deposition of 7-rays for two alternative models for subluminous SNe la: a 
sub-Chandrasekhar edge-lit detonation of a 0.7 Mq C+O WD (model 2 by Livne & 
Arnett 1995), and a pulsating delayed detonation model in a Chandrasekhar C+O WD 
(model WPDl by Woosley, this volume). The deposition profile is calculated at different 
times after the explosion 



^{t) = (6.76 X 10^ D^{t) + 2.72 x 10® Dp{t)) (e"*/^^" - 6"*/^^') 

+ 3.91 X 10^ 1)^(^)6-*/^^' erg g'^ s"^ 

(4) 

where rjvi=8.8 days and rco = 111.26 days are the e-folding times for ra- 
dioactive decay of Ni and Co respectively, and D-y and are the deposition 
functions of 7-rays and e"*", which are decreasing functions of time. Figure 
1 shows D-y for a couple of SNe la models. 

Positrons start to play a major role in the E input as soon as ^ 
0.1. At these low optical depths (achieved around 250 days after explosion 
in a 1 M0 ejecta), most of the 7-ray photons escape. An important issue 
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turns then to be the fate of the e"*" emitted in the ejecta, and the fraction of 
escape and deposition of the kinetic energy of those e"*". The most energetic 
e"*" will succeed escaping the ejecta. The e"*" emitted in the outer layers will 
also have a higher chance to leave the ejecta without becoming thermalized. 
The chance for a e"*" of a given energy to escape or deposit its kinetic energy 
in the ejecta depends very much on the configuration of the magnetic field, 
and requires a numerical evaluation. 

The configuration of the magnetic field of the supernova ejecta is there- 
fore a key subject. How is it possible to ascertain which is the most likely 
configuration of the magnetic field in supernova ejecta? In the first place, 
one can try to track the final configuration from that of the WD. The WD 
can have a dipole magnetic field of primordial origin anchored in the inner 
core, or it can have a dynamo generated magnetic field originated from 
the presence of charges in motion in the envelope of the WD. How is this 
magnetic field preserved after the explosion? It is likely that the convective 
pre-supernova structure will favor a turbulent configuration of the mag- 
netic field, even if it had a dipole configuration in the original WD. From 
the condition that the kinetic energy density of the ejecta should be much 
larger than the energy density of the magnetic field (since the kinetic energy 
density is of the order of the binding energy of the WD), one can estimate 
the intensity of the magnetic field to be: 
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This gives a constraint on the intensity of the magnetic field of B lO^G 
(Chan & Lingenfelter 1993). 

The radial configuration of the magnetic field, as suggested by Colgate et 
al. (1980), would be favored if a dipole magnetic field of the presupernova 
star were combed out by differential expansion of the supernova ejecta, 
and the convective pre-supernova structure had not distorted the assumed 
original dipole configuration. 

A empirical way to evaluate the configuration of the magnetic field is 
to study how much does the late light curve of SNe la depart from the 
full-trapping exponential tail of ^^Co decay. There are reasons to think 
that a radially combed out magnetic field would give too fast a declining 
bolometric luminosity as compared with the observations. 

4. Positrons 

The energy spectrum of the e"*" has the form: 



x{e)^k{Z,e){eo-efeVe^ - I 



(6) 
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where e is the (total) positron energy in nieC^ and 

eo = E-raax/nieC^ + 1 (7) 

k(Z,e)= ) (8) 

with 

Ze^ Za , , 

where Z = 26, f is the speed of the positron, and a is the fine-structure 
constant. 

The positrons lose their energy in the supernova ejecta mainly by ion- 
ization losses. Positrons emitted at a time and place deep in the ejecta 
where the local density is high will likely lose their kinetic energy and ther- 
malize. If a e"*" is emitted at shallow depths or at very late times, such that 
it cannot completely slow down before leaving the ejecta, it will enter the 
interstellar medium as a nonthermal particle. If the density is high enough, 
the thermalized e"*" will tend to annihilate giving two 511 keV photons — 
case the annihilation occurs with a free electron, or they can give either 
two or three 7-ray photons — case the annihilation occurs via positronium. 
The transport and deposition of those 7-rays is usually treated jointly with 
that of the 7-rays emitted in the decay of ^^Co. 

The positron energy loss, per unit column depth, X, due ionization of 
atoms is (Heitler 1954; Blumenthal & Gould 1970; Gould 1972): 

where tq is the classical electron radius, m„ is the atomic mass unit, and 
Z and A are, respectively, the effective nuclear charge and atomic mass of 
the ejecta material. The n(_E) factor is (Berger & Seltzer 1954): 



^{E) = In —, ^ +- n 2- — X — + + -77 + -^^ 

^ ' \ I/m,c^ 12 [2^(7+1)^(7 + 1)2^(7 + 1)2 

(11) 

where 7 = l + E/niec'^ is the Lorentz factor of the relativistic positron, and 
/ is the effective ionization potential for the ambient atoms in the ejecta. 
A semiempirical formula for the ionization potential is (Roy & Reed 1968; 
Segre 1977): 
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I = 9.1z(l + j§^)eV (12) 

Due to the weak dependence of n(E) on /, the above approximate for- 
mula for I is accurate enough for the practical calculation of the energy 
loss. 

Synchrotron losses by the e"*" in the presence of the magnetic field, 
bremsstrahlung losses, and losses due to Compton scattering off photons 
are negligible as compared with those arising from ionization/excitation 
(see Chan & Lingenfelter 1993). Thus, these two last processes are domi- 
nant in determining the distance travelled by the e"*" before losing its kinetic 
energy. 

One can give some typical values, for positrons of different energies, of 
their stopping distance, dg, due to impact ionization and excitation: 

4 ^ - ^ (A) %T + 0-15)-^ cm (13) 

—dL/dX p ymgC^ J Z Iq 



TABLE 1. Typical values for the stopping distance and the 
fraction of the envelope travelled by e"*" of various energies 





1 keV 


10 keV 


100 keV 


1 MeV 


de 


9.5 X 10" 


5.5 X 10" 


3.8 X 10" 


3.0 X 10" 




3.6 X 10"^ 


2.1 X 10"^ 


0.15 


0.11 



de: stopping distance (in cm) 

{ = de/Renv (300 d) 



4.1. TRANSPORT AND ESCAPE OF POSITRONS 

When the positron mean free path is very small as compared with the 
characteristic radius of the supernova ejecta, it is possible to evaluate the 
transport of positrons using the diffusion approximation. This is the case 
when there is a turbulent magnetic field which confines the trayectories of 
the positrons along its lines. The case of the radially combed out magnetic 
field is the opposite: the mean free path of the positron is very large and the 
diffusion approximation is no longer valid: the transport equation has to be 
solved. For the case where the positron has zero diffusion mean free path, 
there is a critical time after explosion, tc(mi,ji), such that positrons born 
at mi with energy at times ti > tc cannot slow down to thermal energies. 
These positrons survive in the eject "fast", nonthermal population. 
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The time tc is defined from the energy loss equation: 

dE 



T{E) 

which is: 



pvdt , (14) 



, p(mi,t' )dt' = —nieC / j^=^= dj (15) 

t, J-,, r(7mec2)^72 - 1 

From the above relation and the homology relation p oc t~'^ , one calculates 
the time, ij, at which a positron, born at time ti, will have lost essentially 
all of its kinetic energy and become "thermalized" (see Chan & Lingenfelter 
1993, for more details): 

tf=i-^--\ , for t, < (16) 
where the critical time after explosion is: 



X / TTT ox / o -, «7 



^'^"""^''^ M \dmj^/\lr r(7mec2)v^^ 

With the decay rate: 



-1/2 



(17) 



R{t) = -(e-*/^^° - e-*/^~0 (18) 

[TCo - TNi) 

the fraction of fast positrons surviving at nii in the ejecta is: 

ffim,) = I de,x{ei)[— exp{-tc{mt,e,)/TNt) (19) 

TNi - TCo Jl TCo 

—exp{-tc{mt,e,)/Tco)] (19) 

x(e) being the positron spectrum (eq. 6). 

In Figure 2 we show the deposition function of the e"*" for the turbulent 
configuration of the magnetic field in different models from 200 to 1000 
days after the explosion. Values for the epoch of interest to our calculations 
are given in Table 2. Chandrasekhar models tend to trap more efficiently 
the energetic e"*" , achieving lower escape fractions than sub-Chandrasekhar 
models. 
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Figure 2. Deposition of energy from e"*" for various models 



TABLE 2. Deposition of e"*" for different models^ 







W7 


DD4 


WPDf 


LAM2 


NIDD 




(300(i) 
(350(i) 


fOO. 
99.9 


99.9 
99.9 


94.9 
94.9 


89.7 
84.6 


79. f 
7f.5 



^Models are: deflagration model W7 of Nomoto, Thielemann & Yokoi (f984); delayed 
detonation model DD4 of Woosley & Weaver (f995); pulsating delayed detonation model 
WPDf of Woosley (this volume); He-detonation of a 0.7 M© WD (model 2 of Livne & 
Arnett f995), and He-detonation of a f .f M© WD (Nomoto f995, here called NIDD) 



When the positrons have essentially infinite diffusion mean free paths and 
they spiral radially through the ejecta, the local fast particle survival frac- 
tion has to be determined numerically. 
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5. Nonthermal Processes 

The Compton scattering of 7-rays has created a distribution of energetic 
e~ called primaries. The energy of the primaries created by Compton scat- 
tering off electrons is given by: 

- COS0) 

" mc2 + E^{1 - cose) 

Those primary e~ are able to give rise to impact ionizations, thus pro- 
ducing a cascade of secondary e~. The typical lower energies of the sec- 
ondaries do not favor further ionization, but they can excite nonthermally 
the ions in the ejecta or lose their energy through Coulomb scattering off 
thermal e~ . 

Treating the energy losses as a continuous process in the case of a plasma 
dominated by Fe (and Fe-like) ions, one can estimate the ionization rate 
due to impact ionization by primaries, secondaries, and the excitation rate. 

The generation of secondaries is given by: 

SiE,^ Es) = \ ,^ (21) 

J ^an-i 1^1 + (f 

where I is the ionization potential and J is an experimentally determined 
parameter varying with the element (Opal et al. 1971), and somewhat lower 
than the ionization potential. 

The cross sections for the different processes competing in the share of 
energy are of the form: 

, , ln(E/I) , , 

cj,,.,{E) a -^Yl ^^^^ 

for ionization and excitation (Lotz 1967), and 

E \lnd 



for Coulomb scattering, where d depends on Ue and Tg. 
Using the available cross sections, one can derive the ionization rate by 
primaries and secondaries. The rates of impact ionization are calculated by 
evaluating the probability of impacts leading to the final result (ionization 
of atom/ion 1) against all other possible paths through which the primary 
(or secondary) can lose their energy. For instance, the impact ionization by 
primaries leading to ions 1 can be evaluated as: 
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^f = n,y:r^dE (24) 

where aik are the corresponding cross sections for the different transitions 
(Axelrod 1980; Swartz 1989, 1991), n,p is the rate of production of primaries 
of energy E,p, and L*{E) is the luminosity loss function which contains 
the sum of the different channels into which the energy of the primaries 
can flow: ionization and excitation in all possible transitions, and Coulomb 
scattering. 

Excitation and ionization rates by secondary e~ is calculated after ob- 
taining the energy distribution of secondaries. The rate at which the energy 
is deposited in heating is evaluated in order to calculate the electron tem- 
perature profile. 



6. The timescales of the radiation processes 

The main atomic timescales involved in the emission at these phases are 
those governing the ionization balance and the thermal balance (since 
excitation-deexcitation processes occur in a very short timescale, steady 
state holds). The recombination timescale is: 

Tree ~ (25) 

where a is the recombination coefficient for a given species and ng is the 
electron density. 

This can be compared with the coUisional excitation timescale, given 

by 

1 , , 

Texc ~ JT- (26) 

Typically, Te^c is much smaller than Tjon, and this ensures the possibility 
of evaluating the statistical equilibrium equations for each ion without the 
need of coupling all energy levels of five or more ionization stages. 
Thermal balance is governed by the dynamical cooling time: 

Tcool = (27) 

a 



and the radiative cooling time: 



Trad = -r 28 
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Ionization balance and thermal balance will be in steady state as long 
as the mentioned timescales will be lower than the evolutionary timescales. 
These are: 

T 56Co III. 2Q days (29) 
the radioactive decay timescale, and the expansion timescale 

n t 

The radioactive timescale reflects the changes in the input energy from 7- 
rays and positrons, and the expansion timescale reflects the changes due 
to homologous expansion and the subsequent decrease in density of the 
ejecta. At around 200-300 days after the explosion, typical number density 
n of the ejecta is 10^-10^ cm and the electron density is of the same 
order. However, around 500 days Ue ~ 10"^ in the central part of the SNIa 
ejecta and r^ec ^ 100 days for Fe ions. Since this timescale is signiflcantly 
longer than the evolutionary timescales governing the changes in the neb- 
ula (radioactive timescale and expansion timescale), then the steady state 
hypothesis for ionization processes breaks down and steady state balance 
is lost. The atomic timescales for recombination processes start to become 
long at about 500 days after explosion. Those for cooling processes remain 
short for quite longer times. 

7. Ionization and rate equations 

In a nonthermal nebula, the ionization and recombination processes occur 
in a longer timescale than excitation and de-excitation, which modify the 
energy level populations in a shorter timescale. Thus, the ionization balance 
can be computed flrst to determine the ion fractions. 

The processes to be included to calculate the ionization balance with 
some degree of accuracy are nonthermal impact ionizations by primaries 
and secondaries, coUisional thermal ionization, and photoionization by emis- 
sion in recombination cascades. 

Since the kinetic energy of the electron is generally much lower than the 
ionization potential (kT/Ip ^ 1), the recombination photons will tend to 
be emitted very near to the photoionization edge. In this case, and if there 
were just a single ionic species, a photon emitted by recombination to the 
ground state of a given ion would be immediately reabsorbed. In the SN 
ejecta the emitted photon can be reabsorbed in other ionizing transitions 
of other species. Terms to account for such interactions can be inserted in 
the same equation of ionization balance. Such equation thus reads: 
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dt 



[rf + rf + Cj{T) + R,k\ n,{t) + a,{T)n,n,+^{t) 



f311 



(in steady state balance for ionization, it turns out to be dni(t)/dt = 0), 
where Ff and Tf are the impact ionization rates by primaries and secon- 
daries and 

R^knt = n, / J^—^dv (32) 

Here, however, the specific intensity of the radiation field Jj, is only due 
to the emission of photons originated in recombination cascade decays and 
the background continuum radiation field is negligible iJ^(^cont) — 0)- 

It is justified to follow the recombination cascade in terms of the possi- 
ble fates of the photons emitted in those transitions. Photons which are not 
going to be absorbed on the spot will likely escape further absorptions: re- 
absorption in permitted transitions has a low optical depth, only enhanced 
near the resonance region, and reabsorption to ionization has a larger prob- 
ability to occur near the blue threshold, where the cross section is larger, 
thus also close to the region where the photon was formed. Here we follow 
an approximate trapping probability approach which links together all ions 
and the possibility of recombination giving photons able to ionize other 
species. It is an escape probability treatment which can account, however, 
for the effect of recombination to states other than the ground state for all 
kinds of ions. 

The total cross section for a photon z/ to be reabsorbed by an ion is 
given by: 

= E ^-^^ (33) 

k=j 

and the optical depth for absorption is: 

Tj = a^^n^r (34) 

Let ^ikj be the probability of a recombination to the a level of ion j 
giving a photon of a frequency z/^, able to give rise to ionization of ion i. 
This parameter is well known for the case of recombination to ground state. 
In general, it is given by: 

^^k, = ^ (35) 



14 



where aj^j is the total cross section for ionization of any ion susceptible 
to be ionized by this photon of frequency z/^, and the function $ enters 
modifying the ionization stage. 
The ionization balance reads: 



{rf + rf + + ^ J2 <i>.fc,«,+i(rKn,+i 

+ (1 - ^tti)a^+i(T)nen^+i 



j i'k>i'i- 



-(1 - ^^-i^^-i^^-i)a,{T)nen, = (36) 

where Vf is the ionization rate due to primaries, Vf is the ionization rate 
due to the secondaries, and C'i (T) is the coUisional ionization rate. 

Such a treatment can be improved by evaluating the flows in the recom- 
bination cascades by independent calculations. Monte Carlo evaluation of 
trapping effects can also help to evaluate the validity of the escape proba- 
bility approach. 

An additional term entering into the ionization balance is the charge 
transfer term. It is not included here due to the unknown rates of the main 
charge transfer reactions involving Fe ions. As a general rule, though, charge 
transfer would tend to modify the fractions of the less abundant ions (the 
intermediate-mass ions) but produce only a small perturbation in the main 
components (Fe-peak ions) (Fransson 1994; Swartz 1994). 

When evaluating the population of the energy levels, a special atten- 
tion is required to build up detailed models for the first five ions of Fe. 
Nonthermal excitation should also be included in the solution of the rate 
equations: in the He I case, due to the energies of levels above the ground 
state, nonthermal excitation can play an important role. 



8. Normal, subluminous, and overluminous SNe la 

The number of data on SNe la studied at these late phases has increased 
over the last years (Kuchner et al. 1995). The observed spectra of SNe la 
show characteristic forbidden lines of Fe"'""'" and Co"'""'" evolving in intensity 
through the late phases. Those observations have shown that the decrease 



15 



of the emission of the lines at 5800 A and 6500 A is consistent with the 
56Qo ^ 56pg (Kuchner et al. 1995). 

There is a diversity among SNe la at late phases. As evidenced by 
the analysis of the spectra, there is a range of velocities over which the 
forbidden line emission of Fe"*" and Fe"'""'" extends, among various super- 
novae. The brighest SNIa, i.e. SN 1991T, has shown the broadest forbidden 
emission both of Fe and Co ions, whereas the faintest ones show narrow 
forbidden emission. A range of ^^Ni masses and electron temperatures has 
also been derived for different SNe la (Ruiz-Lapuente & Filippenko 1993). 
These studies suggest that the amount of ^^Ni synthesized in the explo- 
sions range from 0.4 to 0.8 M©. "Normality" would correspond to about 
0.6 M0, whereas the overluminous SN 1991T would be located at the top 
end (0.8 Mq) and the moderately underluminous SN 1986G at 0.4 Mq 
(Ruiz-Lapuente & Lucy 1992). The very subluminous SN 1991bg can be 
placed at the lower end of the variation range: this SNIa seems to have syn- 
thesized an amount of ^^Ni as low as 0.1 Mq (Ruiz-Lapuente et al. 1993, 
hereafter RJCF93). The spectra of subluminous SNe la indicate differences 
within the class of subluminous SNe la: SN 1986G shows a [Ni II] A 7378 
emission, whereas both SN 1991bg and SN 1991F do not show such emis- 
sion and present a lower ionization stage. SN 1991F can be considered a 
subluminous SNIa of the 91bg-type (Gomez & Lopez 1995). It is possible 
that other subluminous SNe la like SN 1992bc, SN 1992K, and SN 1971J 
would have spectra similar to SN 1991bg. Such observations were never 
performed, however. Thus, the late-time phases of subluminous SNe la are 
poorly documented as compared with those of "normal" SNe la, such as 
SN 1990N, or SN 1992A. 

Limits on the fraction of Fe° and Fe"*""^ have been established respectively 
from the infrared (Spyromilio et al. 1994) and the UV (Ruiz-Lapuente et 
al. 1995, hereafter RKPC95) emissions. By studying the nebular emission 
of SNe la it can be concluded that the dominant ionization stage is Fe"'""'". 
The ionization stage is lower in the most subluminous SNe la. 

The presence of intermediate-mass elements in the Fe-peak region of 
the ejecta of the faint SN 1991bg suggested that the burning took place 
at lower densities than in "normal SNe la". Detonation models of sub- 
Chandrasekhar WDs gave a resonable account of the observed spectrum 
(RJCF93). We will also discuss the nebular spectra of pulsating delayed 
detonation models as alternative candidates to explain subluminous SNe la 
(§10). 
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9. Diagnostics of the supernova density and mass 

Several forbidden transitions of Fe"*" and Fe provide interesting diagnos- 
tics of the electron density ng and electron temperature Tg of the supernova 
ejecta. We would like to outline here how this can be useful for the present 
debate between sub-Chandrasekhar and Chandrasekhar models. The emis- 
sivities of the lines in a variety of astrophysical situations clearly indicate 
the physical conditions of the region in which they originate. The approach 
of infering the electron density profile rae(r) of the ejecta through the emis- 
sivity profile, and extracting reliable information on n(r) is feasable and can 
carry little error if it is done using the whole information on line emissivities 
together with the most updated results of the Opacity Project (Pradhan & 
Berrington 1993; Bautista & Pradhan 1994). 

There is a number of emission lines which can be used in the density 
diagnostics. As in other nebulae, the population of energy levels is going 
to evolve as the supernova density decreases due to expansion. The upper 
levels become underpopulated (with departure factors of 10~^-10~"^). For 
some models, this produces some characteristic transitions to show much 
fainter emission than observed for a given phase. A clear diagnostics of 
low rie is associated to the faintness of forbidden emission at A 5200 A, 
A4300 A and A5000 A, due to Fe+ a^F-b^P and a^F-a^H multiplets. The 
lower energy terms of these transitions can be significantly depopulated 
if rie is low. Interesting line ratios are A5262/A8617 and those involving 
the emission at 4400 A ([Fe II] A 4416 and other transitions) as compared 
with the emission in the red and the infrared. Those transitions can help 
to discriminate the central density of a given supernova. Therefore, this 
can help to track down the mass. Here we show the evolution of some of 
these ratios for models with different masses: in Figure 3 we illustrate how 
the ratio A5262 to A8617 evolves with time, for different explosion models 
involving different WD masses, and for two phases: 270 and 300 days after 
explosion. 

There are other transitions of interest for which collision strength cal- 
culations are available. The infrared transitions help as indicators of the 
electron density, in particular the [Fe II]A 8617 to [Fe II]A 1.257 //m line 
and, in general, the ratios of lines arising from different energy levels. This 
approach is very useful, since it overcomes other uncertainties, thus re- 
ducing the source of error to the accuracy of collision strengths for the 
transitions considered. The total mass of the ejecta can be tracked down 
accurately from the electron density profile, due to the univocal relation- 
ship between mass of the exploding WD and central density at any time 
after explosion. 
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Figure 3. Ratio of two characteristic lines for density diagnostics 



10. Models 

Sub-Chandrasekhar models have a number of advantages when the vari- 
ation of properties among SNe la is to be explained. They can account 
easily for the variation of luminosities and rates of decline of the light 
curves (Phillips 1993; Riess, Press & Kirshner 1994; Livne & Arnett 1995; 
Woosley & Weaver 1994). They also explain in an easy way the correlation 
of these properties (rate of decline of the light curve and maximum bright- 
ness) with galaxy type, favoring fainter SNe la in early type galaxies as a 
result of the predominance of older and less massive white dwarfs among 
the exploding objects (Ruiz-Lapuente, Burkert & Canal 1995). 
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SN 1994D in NGC 4526 
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Figure 4- Comparison of the spectrum of a sub-Chandrasekhar model (model 6 by 
Livne & Arnett 1995) with the one of the Chandrasekhar deflagration W7 (by Nomoto, 
Thielemann & Yokoi 1984) 



Does, however, the range of sub-Chandrasekhar type of explosions cor- 
respond to what is really observed? One would expect, if sub-Chandrasekhar 
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Figure 5. Electron density profiles for different models. There is a significant range of 
central electron densities in models of exploding WDs of various masses at a given phase. 
This property can help us to discriminate between models 



mass WDs exploded in a continuous range from 0.6 to 1.0 M©, to have a 
variety of late-time spectra which has not been seen from the observations 
so far. Here we judge on the basis of the assessment of ID He-detonation 
models (Nomoto 1995, private communication; Woosley & Weaver 1994, 
and our own calculations) as well as 2D calculations (Livne & Arnett 1995). 
Both ID and 2D calculations seem to give similar results. 

The central and upper range in mass of the sub-Chandrasekhar C + 
explosions give expanding ejecta of density too low as compared with that 
inferred from observations. The signature of such a low density is the ab- 
sence of some characteristic emission arising in upper energy levels of Fe"*" , 
since those levels become less populated by collisions in lower density ejecta. 
What happens with the lower mass end of this set of models? The observed 
spectrum of SN 1991bg and the rate of decline of its light curve seems 
to be reasonably well explained by the sub-Chandrasekhar models. Both 
ID and 2D He detonation models do, however, yield a bit too much ^^Ni 
as compared with what is found at late phases in this supernova. The al- 
ternative explanation to this supernova, within the Chandrasekhar-mass 
WD paradigm, can be supplied by pulsating delayed detonation models 
(see Hoflich, Khokhlov & Wheeler 1995, and Woosley, this volume). Model 
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Figure 6. A possible model for the subluminous SN la SN 1991bg. The solid line shows 
the model spectrum, and the dashed line the observed one 

WPDl by Woosley (this volume) would be the closest one, among Chan- 
drasekhar models, to SN 1991bg. Ways to discriminate between the sub- 
Chandrasekhar and Chandrasekhar possibilities are: (a) through the spec- 
tral evolution at late phases, i.e. faster in the low WD model than in the 
pulsating delayed detonation; (b) through the ionization stage shown at 
these phases, with [Fe III] emission in the low WD explosion model and 
dominant [Fe II] in the faint WPDl model; (c) through the presence or 
not of the [Fe II]A5557 emission (and other emissions associated to mul- 
tiplets a'*D-b^P and a'*P-b'*D), which become prominent if is about 
10'' cm~"^ and Fe"*" is the dominant ionization stage, conditions found in 
Chandrasekhar-mass explosions which synthesized a low amount of ^^Ni 
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(such as some pulsating delayed detonations). Such emission, in contrast, 
is not expected to be important in a sub-Chandrasekhar envelope, due to 
the lower ng. Both types of explosions have in common producing faint 
red SNe la at maximum, which are also redder and fainter at late phases 
than "normal" SNe la. A more complete set of observations is needed to 
ascertain whether 91bg-like events correspond to the lower edge of sub- 
Chandrasekhar models or to the fainter set of Chandrasekhar pulsating 
delayed detonations, having in view the above discriminating clues. 

11. Late light curves 

As it has long been argued, the light curves at late phases constitute one 
more element in the discussion of the total ejected mass in Type la SNe 
(Colgate 1991). At the phase where the thermalization timescale of the 
radioactive input is short and work to accelerate the envelope is negligible, 
the luminosity describes the convolution of the exponential decay of ^^Co 
with the increasing transparency (and escape) of both 7-rays and e"*". 

From observations, it is seen that the monochromatic light curves of 
SNe la follow a steeper decline than the full trapping curve of e"*" at the 
times where these dominate the late luminosity (and 7-rays fully escape). 
Very few observations allow to reconstruct bolometric light curves: the re- 
constructed bolometric light curve of SN 1992A (Suntzeff 1995) shows a 
slight departure below the full trapping curve, whereas in the faint SN la 
SN 1991bg the departure from the full trapping curve (Turatto et al. 1995) 
seems larger. 

Given the nonthermal, NLTE character of the radiating SN ejecta, and 
the timescales involved in the different processes, the derivation of the 
monochromatic light curves has to be done through the solution of the spec- 
tral emission at all wavelength ranges. Here, we show monochromatic light 
curves up to the phase where time-dependence in ionization and thermal 
balance starts to play a role. To sample the light curve we have performed 
spectral calculations every 15-20 days in the tail of the light curve (the 
colors are given in the B, V, R Johnson system). 

Chandrasekhar models for the turbulent configuration of the magnetic 
field show almost full trapping of the energy of e"*" up to 1000 days, whereas 
edge-lit detonations (sub-Chandrasekhar explosions) predict a departure 
of the bolometric luminosity below the full trapping line already at 200 
days, for the same magnetic field configuration. In the case of a radially 
combed out magnetic field, both Chandrasekhar and sub-Chandrasekhar 
models experience an earlier departure from the radioactive tail, which is 
likely to be too large as compared with most SNe la ("normal" ones). Figure 
7 shows the result of making the comparison of monochromatic theoretical 
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light curves with observed ones. It is far more illustrative than the compar- 
ison of the bolometric light curves, since the monochromatic decays of the 
light curves inform on how the luminosity in the different spectral bands is 
distributed and this is much more directly linked to the "goodness" of the 
underlying model. The agreement of the model in velocity, composition and 
density space can be, however, only ascertained through the full spectral 
comparison. 




Figure 7. The solid lines show the late B V R light curves of model WPDl (Woosley, 
this volume), the dashed lines display the same color light curves for model 2 of Livne 
& Arnett (1995). The solid points correspond to the observations of SN 1991bg (Turatto 
et al. 1995). The doted line in the central panel is the visual light curve of SN 1992A. 



12. Final remarks 

There is a large number of diagnostics that we can use to obtain infor- 
mation on the validity of a given SN model. Late-time spectra and light 
curves provide an ideal frame to discuss the mass of the ejecta, kinematics, 
and nucleosynthesis of thermonuclear SNe. Through the exploration of the 
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ratios of different lines and tlieir dependence on tlie electron density, we 
gain insight on whether SNe la are sub-Chandrasekhar or Chandrasekhar- 
mass explosions. The "present" set of sub-Chandrasekhar explosions seems 
to give too low density profiles: about 300 days after explosion the emission 
coming from upper energy levels is significantly depleted as compared both 
with more massive models and with the observations. In Chandrasekhar- 
mass models, emission from excited levels (such as the emission at 4400 A) 
is well reproduced. Red faint SNe la at late phases are expected both from 
low-mass sub-Chandrasekhar explosions and from Chandrasekhar pulsating 
delayed detonations with little (^0.1 Mq) ^^Ni. Electron density diagnos- 
tics can help to discriminate between these two possibilities. 

A related way to pin down the correct thermonuclear supernova model 
is through the study of the light curve tails. Chandrasekhar models trap sig- 
nificantly the 6+ and will give bolometric light curve declines close to the full 
trapping line drawn by the exponential decay of ^^Co. The bolometric light 
curve of sub-Chandrasekhar models tends to fall below the full-trapping 
line of ^^Co decay. Monochromatic light curves (B, V, & R) can experi- 
ence steep declines in some SNe la models, as for instance in model 2 of 
Livne & Arnett (1995), or in the ^^Ni-poor pulsating delayed detonations 
(such as WPDl of Woosley, this volume). If dumpiness is present in the 
supernova ejecta at this phase, escape will be enhanced and the outcome 
can also be light curves declining faster at late phases. In pulsating delayed 
detonation models of weak explosions, the formation of dust grains (due 
to lower temperature), can lead to additional departures in the monochro- 
matic light curves from those for full-trapping. If there is mixing between 
the region where intermediate-mass elements are present together with Fe 
(from the previous ^^Ni ) with the region, the formation of molecules 
such as MgSiOs or Fe304 would be favored in models such as WPDl (the 
electron temperature is below 2000 K at about 300 days after explosion in 
those regions). 

A more systematic observational coverage of late-time spectra of ther- 
monuclear SNe la should be done. So far, the sample of "normal" SNe la 
has been well tracked into the late phases. However, late-time observations 
of the subluminous SNe la have been very unfrequent. It is important to 
realize that a good nebular tracking of those supernovae will help to dis- 
criminate whether WDs below the Chandrasekhar mass do explode or not. 

I would like to thank especially Eli Livne and Stan Woosley for providing 
me with their hydrodynamic models, and for related discussions. My thanks 
go as well to Leon B. Lucy and Anil Pradhan for informations on available 
atomic data and discussions on radiation processes in SN nebulae. 



24 



References 

1. Axelrod, T.S. 1980. PhD thesis. U of California 

2. Bautista, M. A. & Pradhan, A. K. 1995, A & A (in press) 

3. Belger, M.J. & Seltzer, S.M. 1964, Tables of Energy Losses & Ranges of Electrons & 
Positrons (Washington, DC: NASA) 

4. Blumenthal, G.B., & Gould, R.J. 1970, Rev. Mod. Phys., 42, 237 

5. Chan, K. & Lingenfelter, R.E. 1993, ApJ 405, 614 

6. Colgate, S.A. 1970, ApSS 8, 457 

7. Colgate, S.A., Petschek, A.G., & Kriese, J.T. 1980, ApJ, 237, L81 

8. Colgate, S.A. 1991, in Supernovae, ed. S. E. Woosley (Springer- Verlag), 585 

9. Eastman, R.G. & Pinto, P.A. 1993, ApJ 412, 731 

10. Fransson, C. 1994, in Proceedings of the Les Houches on Supernovae, ed. S.A. 
Bludman, R. Mochkovitch & J. Zinn-Justin. 

11. Fransson, C. & Chevalier, R.A. 1989, ApJ 343, 323 

12. Gomez, G. & Lopez, R. 1995, AJ 109, 737 

13. Gould, R.J. 1972, Physica, 60, 145 

14. Heitler, W. 1954 The Quantum Theory of Radiation, (3d ed.; Oxford: Clarendon 
Press) 

15. HSflich, P., Khokhlov, A.M., & Wheeler, J.C. 1995, ApJ 444, 831 

16. Kuchner, M. J., Kirshner, R.P., Pinto, P.A. & Leibundgut, B. 1994, ApJ 426, L89 

17. Livne, E. & Arnett, D. 1995, ApJ 452, 62 

18. Lotz, W. 1967, ApJ Suppl 14, 207 

19. Meyerott, R.E. 1978, ApJ 221, 975 

20. Meyerott, R.E. 1980, ApJ 239, 257 

21. Nomoto, K., Thielemann, F.-K. & Yokoi, K. 1984, ApJ 286, 644 

22. Opal, C.B., Peterson, W.K., & Beaty, E.G. 1971. J. of Chemical Phys. 55, 4100 

23. Phillips, M.M. 1993, ApJ, 413, L105 

24. Pradhan, A.K. & Berrington, K.A. 1993, J. Phys. B, 26,157 

25. Riess, A. G., Press, W.H., Kirshner, R.P. 1995, ApJ 438, L17 

26. Roy, R. R. & Reed, R.D. 1968, Interactions of Photons and Leptons with Matter 
(New York: Academic) 

27. Ruiz-Lapuente, P. & Filippenko, A.V. 1993, in Origin and Evolution of the Ele- 
ments, ed. N. Prantzos, E. Vangioni-Flam, & M. Casse (Cambridge: Cambridge Univ. 
Press), 318 

28. Ruiz-Lapuente, P. & Lucy, L. B. 1992, ApJ, 400, 127 

29. Ruiz-Lapuente, P., Jeffery, D. J., Challis, P. M., Filippenko, A.V., Kirshner, R.P., 
Ho, L.C., Schmidt, B.P., Sanchez, F. & Canal, R. 1993, Nature, 365, 728 

30. Ruiz-Lapuente, P., Kirshner, R.P., Phillips, M.M., Challis, P.M., Schmidt, B. P., 
Filippenko, A.V. & Wheeler, J.C. 1995, ApJ 439, 60 

31. Ruiz-Lapuente, P., Burkert, A. & Canal, R. 1995, ApJ, 447, L72 

32. Segre, E. 1977, Nuclei and Particles (New York: Benjamin / Cummings) 

33. Spyromilio, J., Meikle, W.P.S., Allen, D.A., & Graham, J.R. 1992. MNRAS 258, 
53p 

34. Suntzeff, N. B. 1995., in lAU CoUoq. 145, Supernovae and Supernova Remnants, 
ed. R. McCray (Cambridge: Cambridge Univ. Press), in press 

35. Sutherland, P. G. & Wheeler, J.C. 1984. ApJ 280, 282 

36. Swartz, D.A. 1989. Ph. D. Thesis. University of Texas, Austin 

37. Swartz, D.A. 1991. ApJ 373, 604 

38. Swartz, D.A. 1994. ApJ 428, 267 

39. Swartz, D.A., Sutherland, P.G. & Harkness, R.P. 1995. ApJ 446, 766 

40. Turatto, M. et al. 1995, MNRAS, submitted 

41. Woosley, S.E. & Weaver, T.A. 1994, ApJ 423, 371 

42. Woosley, S.E. & Weaver, T.A. 1995, in Proc. Les Houches, Session LIV, Supernovae, 
ed. S.A. Bludman, R. Mochkovitch, & J. Zinn-Justin (Amsterdam: Elsevier), 63 



